Cosmic Microwave Background Polarization Signals from Tangled Magnetic Fields 
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Tangled, primordial cosmic magnetic fields create small rotational velocity perturbations on the 
last scattering surface (LSS) of the cosmic microwave background radiation (CMBR). For fields 
which redshift to a present value of Bo = 3 x 10 -9 Gauss, these vector modes are shown to generate 
polarization anisotropies of order 0.1/iif — 4fiK on small angular scales (500 < I < 2000), assuming 
delta function or a power law spectra with n = — 1. About 200 times larger signals result for n — 2 
spectra. Unlike inflation generated, scalar modes, these signals are dominated by the odd parity, 
B-type polarization, which could help in their detection. 
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Magnetic fields in astronomical objects, like galaxies, 
could grow by the amplification of small seed fields by 
turbulent dynamo action Qj. However, the need to pro- 
duce magnetic helicity in galaxies seems to severely con- 
strain the efficiency of such dynamo action ||. Alter- 
natively, the galactic field could be of primordial origin 
||, although, there is no compelling mechanism for pro- 
ducing the required field ^ . A primordial field that red- 
shifted to a present value of ~ 10 -9 Gauss, tangled on 
galactic scales, could also significantly affect galaxy for- 
mation It is of considerable interest, therefore, to 
find ways of constraining or detecting such fields 
Observations of anisotropies in the CMBR, provide a po- 
tentially powerful constraint. Indeed, the CMBR tem- 
perature isotropy can be used to place limits, of order 
several nano-Gauss, on both the uniform || and tangled 
components of the magnetic fields |]|n|. However, tem- 
perature anisotropics in inflationary models of structure 
formation, will be dominated by "non-magnetic", scalar 
modes. So the detection of a magnetic field induced 
signal is likely to be difficult, except possibly on scales 
smaller than the Silk damping scale P). Here we point 
out the advantage of using alternatively, the polarization 
anisotropy. Note that scalar perturbations only produce, 
what is known as E-type polarization anisotropy. How- 
ever, as we show here, tangled magnetic fields which drive 
significant vector perturbations, also lead to a distinc- 
tive, significant and potentially detectable B-type polar- 
ization anisotropy of the CMBR. This could help in sep- 
arating their contribution from scalar contributions, to 
detect/constrain such tangled, primordial fields. 

Polarization of the CMBR arises from the Thomson 
scattering of radiation from free electrons, and is sourced 
by the quadrupole component of the CMBR anisotropy. 
The evolution equations of temperature and polariza- 
tion anisotropy for vector perturbations have been de- 
rived in detail in Ref. in the total angular momen- 
tum representation. We use their results extensively be- 
low. The anisotropy in the temperature and polariza- 
tion is expanded in terms of tensor spherical harmon- 
ics. This enables one to write evolution equations, for 



the moments, o[ m ' ) , £7, ' and B\" L ' , of the temperature 
anisotropy (AT/T), the electric (E-) type and the odd (56) of |l_] 



(m) 



parity, magnetic (B-) type polarization anisotropies, re- 
spectively. Here I stands for the multipole number and 
m = 0,±1,±2, respectively, for scalar, vector and ten- 
sor perturbations. For vector perturbations (m = ±1), 



the magnetic ty 



De contribution dominates the polariza- 



tion anisotropy | 11J|. Its evolution is given by (Eqn. (77), 
(62) and (18) ofp) ), 
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yp — l)(l + 2)ji(x)/2x, with ji(x) the spherical Bessel 
function of order I. The 'visibility function', <?(to,t) = 
k(t) exp[— k{r')dT'\, determines the probability that 
a photon reaches us at the conformal time To if it was last 
scattered at the epoch r. Here k(r) = n e (r)<Jxa(T), n e 
the electron number density, <jt the Thomson cross sec- 
tion, and a(r) the cosmic scale factor normalised to unity 
at the present. We assume a flat universe throughout. 

For standard recombination, g is sharply peaked about 
the time of recombination. We therefore need to calcu- 
late and hence the quadrupole anisotropies around 
this epoch of last scattering. These can be analyti- 
cally estimated using the tight-coupling approximation, 
k/k = kL 1 <C 1. Here L 1 {t) = (k) _1 is the co-moving, 
photon mean free path. First, to leading order in this 
approximation, we have zero quadrupoles, and a dipole 



where is the magnitude of the (vector 



u l — U B ' 

or rotational component of) baryon fluid velocity field, in 
Fourier space. However, to the next order the quadrupole 
is not zero. It is generated from the dipole at the 'last 
but one' scattering of the CMBR. Using the moments of 
the Boltzmann equations for the temperature and polar- 
ization anisotropies (Eq. (60), (63) and (64) of p|), one 



gets 6^ = 
P« = 0^/4 = kL^i 
gives an estimate of B 
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/(3\/3). Using this in Eq. (g) 
and the angular power spectra 
y. We us 
and get 



Cf B due to B-type polarization anisotropy. We use Eq 



to relate Cf B interms of .B, 
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Here we have included an extra factor of 2, since we have 
to sum over the power in both m = + 1 and m = — 1 con- 
tributions. The above expression for C BB is very closely 
related to Eq. (1) of Ref. || (henceforth, Paper I), for 
the temperature power spectrum C; , due to tangled mag- 
netic fields. One can make the same approximations as 
made there, to obtain an analytic estimate of C BB . 

Firstly, it suffices to approximate the visibility func- 
tion as a Gaussian: <?(t ,t) = (27rcr 2 ) -1 / 2 exp[— (r — 
t*) 2 / (2a 2 )], where t* is the conformal epoch of "last 
scattering" and a measures the width of the LSS. Us- 
ing the expressions given in Ref. |l2|, we estimate 7% ~ 
176.2h~ 1 Mpc and a = llM^Mpc. (h is the Hub- 
ble constant in units of 100 km s _1 Mpc -1 . We take 
h = 0.75 henceforth and adopt a baryon density param- 
eter = 0.02ft, -2 .) The dominant contributions to the 
integral over r in Eq. (^) then comes from a range a 
around the epoch r = t*. Further ji(k(ro — t)) picks out 
(k, t) values in the integrand which have fe(ro — r) ~ I. 

Let us also assume that kL^v^ (k, r) varies slowly with 
r (around r ~ r*), and slowly with k around k ~ Z/i?*, 
the regions which contribute dominantly to the integrals 
in Eq. (^). (Here i?* = tq — r»). Then following exactly 
the arguments detailed in paper I, we get for ka << 1, 
the analytical estimate, 
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Here, A 2 = k 3 < \v£\k,T*)\ 2 + {v^ (k,n)\ 2 > /(2tt 2 ) 
is the power per unit logarithmic interval of k, residing 
in the net rotational velocity perturbation. And in the 
other limit, ka >> 1, we get 
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For small wavelengths, Cf B is suppressed by a 1/ka fac- 
tor due to the finite thickness of the LSS. Note that 
in both cases, the polarization anisotropy, AT BB (l) w 
(fcL 7 (r*)/3) x AT(l), where, AT(l) is the temperature 
anisotropy computed in Paper I. 

To evaluate C, BB , one needs to estimate 



>B > 



for a 

general spectrum of magnetic inhomogeneities. We as- 
sume the magnetic field to be initially a Gaussian ran- 
dom field. On galactic scales and above, the induced 
velocity is generally so small, it does not lead to any ap- 
preciable distortion of the initial field ||. So, to a very 
good approximation, the magnetic field simply redshifts 
as, B(x, t) = bo(x)/a 2 . The Lorentz force associated 
with the tangled field is then = (V x b ) x b /(47ra 5 ), 



which pushes the fluid, creating rotational velocity per- 
turbations. Further, the magnetic stresses, say n^, can 
lead to metric perturbations. We focus on scales larger 
than the photon mean-free-path at decoupling, and de- 
scribe the viscous effect due to photons, in the diffusion 
approximation. The Fourier transform of the linearised 
Euler equation for , is given by |i~l"|,p|], 
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Here V(k, t) is the vector component of the metric per- 
turbation (t is comoving proper time), p 7 the photon 
density, pi, the baryon density, and rj = (4/15)p 7 / 7 the 
shear viscosity coefficient associated with the damping 



due to photons, whose mean-free-path is 



L 7 a(t). 



B 



= 1 component (defined as in ref. [fill) 
of PijFj, the rotational part of the Lorentz force. We 
have defined the Fourier transforms of the magnetic field 
as, bo(x) = b(k) exp(ik.x). Since the Lorentz force 
is non-linear in bo(x), this leads to the mode-coupling 
term F(k) = £ p [b(k + p).b*(p)]p - [k.b*(p)]b(k + p). 
The projection tensor, -Py (k) = [£y — kikj/k 2 ] projects 
F onto its transverse components (perpendicular to k ). 

The comoving Silk scale at recombination, Lg — kg ~ 
(Z 7 (i*)t*) 1 / 2 /a(i*) ~ 6.8Mpc, separates scales for which 
the damping term in (|5|), is important (kLs >> 1) from 
those for which it is negligible (kLs << 1) Jjj§- We can 
solve Eq.(|5|) analytically, in these two limits. For kL s < 

1, and when the fluid starts from rest (vg\n) = 0), the 
damping due to the photon viscosity can be neglected 
compared to the Lorentz force. Integrating Eq.(|3|) gives 



value of p 7 , and 5* = (3p b /4p 7 )(r„) - 0.73(fi fc /0.02/i- 2 ). 
The metric perturbation term V is also smaller than 



Gzj(t- t0/(1 + where we have defined 
Gn' = 3Fg^/[167rpo], with Po the redshifted present day 



the Lorentz force driven contribution to , for large 
\k=i/R,- (4) I by a factor ~ (l/24)~ 2 h~ 2 (see p5|); and so makes 
a negligible contribution to C BB , for the small angular 



scales (I > 400) considered here. In the other limit, with 
kL s >> 1, we can use the terminal- velocity approxima- 
tion, neglecting the inertial terms in the Euler equation, 
and balance the Lorentz force by friction. This gives 
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(G^ /k)(kL 7 /5)~\ independent of V. 



We also need to specify the spectrum of the tangled 
magnetic field, M(k). We define, < £>j(k)6j-(q) >= 
<5k,q-Fjj(k)M(/c), where 5k, q is the Kronecker delta which 
is non-zero only for k = q. This gives < b 2 , >= 
2j(dk/k)A 2 (k), where A 2 (k) = k 3 M (k) / (2n 2 ) is the 
power per logarithmic interval in k space residing in mag- 
netic tangles, and we replace the summation over k space 
by an integration. The ensemble average < |w^| 2 >, 
and hence the Cp B s, can be computed in terms of the 
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magnetic spectrum M(k). It is convenient to define a di- 
mensionless spectrum, h(k) = A 2 (fc)/(.Bo/2), where B 
is a fiducial constant magnetic field. The Alfven velocity, 
Va, for this fiducial field is, 



V A 



Bo 



3.8 x 10 _4 B_ 
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where £>_g = (£?o/10 _9 Gauss). Also, as a measure 
of the B-type CMBR polarization anisotropy induced 
by the tangled magnetic field, we define the quantity 
~~(l) = [l{l + l)C BB /2^] 1 l 2 T . 
Note that when ka < 1 (or I < 500), one also 
generally has kL s < 1 0. The resulting AT BB 
can be estimated using Eq.(^|), and ignoring viscous 
damping. A lengthy calculation gives, for such scales, 
AT BB (l) = T (7r/32) 1 / 2 /(fc)[fc 2 ViT,L 7 (rO/3(l + S m )], 
where, I = fci?*. For scales with kLs > 1 (I > 1150), 



AT BB 



we can use Eq.(H), and v 
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similar calculation to that above gives, AT BB (l) — 
T (57r 1 /4/12V2)(/(/c)y|/(fc(7) 1 / 2 ). The function I 2 (k) is 
a dimensionless mode-coupling integral given by H 
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where |(k + q)| = (k 2 + q 2 + 2kqp) 1 / 2 . Putting in nu- 
merical values we estimate for I < 500 and I > 1150 
respectively, 
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Further, for intermediate scales, 500 < I < 1150, an es- 
timate using Eq.(|j), but ignoring viscous damping, gives 
AT BB (l) ~ 0.5^(//800) 3 / 2 /(Z/i?»)(£- 9 /3) 2 . 

If the magnetic spectrum has a single scale, with 
h(k) = kSo(k — fc ), where S D (x) is the Dirac delta 
function, < >— Bq and the mode-coupling integral 
can be evaluated exactly. We find I(k) = (k/ko)[l — 
(fc/2fc ) 2 ] 1/2 , for k < 2k , and I(k) = for larger k. 
So I(k) ~ 1 for k ~ fco, with a maximum /(\/2fco) = 
1. For B_ 9 - 3, one then expects a RMS ATp 33 ~ 
0.2fiK — 1/xiiT, depending on fco and i. We can also 
consider power law spectra, M(k) — Ak n cut-off at 
fc = fc c , where fc c is the Alfven- wave damping length-scale 
|p|,^6[. We fix A by demanding that the field smoothed 
over a "galactic" scale fcc = 1/iMpc -1 , (using a sharp 
fc-space filter) is Bq, giving h(k) = (n + 3)(k/kc) 3+n 
(n > —3). We then find for fc << fc c (as is relevant 
for I < 2000), and n > -3/2, 7 2 (fc) = (28/15)((n + 
3) 2 /(3 + 2n))(fc/fc G ) 3 (fc c /fc G ) 3 + 2 ™. Using this in Eq.(|), 
we find ATjP - 0.16^(l/400) 7 / 2 for I < 500 and 
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FIG. 1. ATp B versus I plot for different model P(k) tak- 
ing B-g — 3. The dashed-dotted line is for the delta function 
spectrum with fco = 0.075 Mpc -1 . The solid curve, is for the 
power law spectrum, with n = — 1, while the dashed curve, 
gives ATp S /I00 for n = 2. For comparison, we also show 
(dash-triple dotted line), the tensor contribution to B-type 
polarization, in a standard cold dark matter model, computed 
with the default parameters of CMBFAST, except changing 
the scalar spectral index to ns = 0.9. The dotted line, gives 
the contribution of galactic foreground emission. 



AT BB - 5^AT(Z/1500) for I > 1150, for n = -1 and 
£?_g ~ 3. Much larger signals result for larger n = 2, 
say, and same B , with AT BB ~ 9.6fiK(l/400) 7 ^ 2 for 
I < 500 and ~ 308.5/iif(Z/1500) for I > 1150. 

To complement the analytics, we have also computed 
for the above spectra, by evaluating the r and 
fc integrals in Eq.(||) numerically (but using analytical 

approximations to v^). The results are shown in Fig- 
ure 1. We see that for Bq ~ 3 x 10 _9 G, one expects 
a RMS B-type CMBR polarization anisotropy of order 
0A/j,K - AfiK for 500 < Z < 2000, for delta function 
or a power law spectra with n = — 1, and ~ 200 times 
larger signals for n = 2 spectra. These values compare 
reasonably with the analytical estimate from Eq. (||) (see 
also [0). Further, both the tensor contribution of typi- 
cal inflationary models ( from CMBFAST 0), and the 
foreground contribution (dominated by galactic dust po- 
larized emission for high frequencies cf. |l9|) are sub- 
dominant to the magnetic field induced signal, for large 
I > 400. Clearly with the sub-micro Kelvin sensitivities 
expected from Planck |2(J, these signals can be detected. 

Earlier work also emphasised the possibility of Fara- 
day rotation and the depolarization of the CMBR due 
to differential Faraday rotation in a tangled magnetic 
field . Note that the average Faraday rotation 
(in radians) between Thomson scatterings is given by 
F = 3B /(27rei/g) w 0.23(B /3 x lQ- 9 G){v ) 2,0GHz)" 2 , 
where v§ is the observed frequency. The CMB could then 
become significantly de-polarised due to this effect, for 
vq < 16AGHz(B Q /3 x lO^G) 1 / 2 , but Faraday rotation 
effects will be negligible for say vq > 40GHz, or the 
higher frequency instruments of the Planck Satellite. 

Note that scalar modes can also be induced by tangled 
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magnetic fields and generate a purely E-type polariza- 
tion. These are however of smaller amplitude than vec- 
tor modes, due to the larger restoring force contributed 
by the radiation-baryon fluid pressure f| . They are also 
strongly damped on scales smaller than the Silk scale, 
Ls, while the vector mode, survives on much smaller 
scales > VaLs @, or larger I (see Figure 1). The ten- 
sor mode can also contribute to AT BB , but one can show 
that their effect is only important at large angular scales. 
Further, AT BB could also be generated, after the uni- 
verse gets re-ionized, though from the recent detection of 
multiple doppler peaks, the optical depth for scattering 
at re- ionization is probably very small p^ |. A detailed 
computation of these effects will be presented elsewhere. 

In conclusion, we have identified here, a new physical 
effect of tangled magnetic fields; that they can produce 
distinctive and potentially detectable B-type polarization 
anisotropy on arc minute scales. From Eq.(|J) and Fig- 
ure 1, we see that a tangled field with B ~ 3 x 10 _9 G, 
induces a RMS B-type CMBR polarization anisotropy 
of order O.lfiK — 900^i£~ or larger, depending on M{k) 
and I. The anisotropy in hot/cold spots could be sev- 
eral times larger, because the non-linear dependence of 
Cf B on M(k) will imply a non-Gaussian statistics for 
the anisotropies (see Paper I). Further in standard mod- 
els all the Cis have a sharp cut-off for I > R*/Ls, due to 
Silk damping but strong damping of Alfvenic modes is 
expected only on scales smaller than VaLs @- Finally, 
since tangled magnetic fields produce predominantly B- 
type polarization, which are also dominant at large I, 
they can be distinguished from those produced by infla- 
tionary scalar and tensor perturbations. Satellite bourne 
experiments like Planck p0[ , with the sub-micro Kelvin 
sensitivities should be able to detect and isolate the ef- 
fects of magnetic fields, using CMBR polarization, if such 
fields indeed play a role in structure formation. 

As this letter was prepared for (re)submission, 
preprints [p3| appeared which have some overlap with 
our work. We thank H. M. Antia and S. Sethi for help. 
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